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We provide a comparison of accurate full potential band calculations of the Fermi surfaces 
areas and masses of MgB2 and ZrE>2 with the de Haas- van Alphen date of Yelland et al. and 
Tanaka et al. respectively. The discrepancies in areas in MgB2 can be removed by a shift 
of a bands downward with respect to ir bands by 0.24 eV. Comparison of effective masses 
lead to orbit averaged electron-phonon coupling constants A CT =1.3 (both orbits), A^=0.5. 
The required band shifts, which we interpret as an exchange attraction for a states beyond 
local density band theory, reduces the number of holes from 0.15 to 0.11 holes per cell. This 
makes the occurrence of superconductivity in MgB2 a somewhat closer call than previously 
recognized, and increases the likelihood that additional holes can lead to an increased T c . 



I. INTRODUCTION 

The discovery of superconductivity in MgB2 near 
40 K by Akimitsu's group jl| and the subsequent 
intense experimental investigation of its properties 
and theoretical exposition of the underlying causes, 
has made it clear that MgB2 is the first member 
of a new class of superconductors. Although inter- 
metallic like the best superconductors before 1986, 
it does not have d electrons, it does not have a 
high density of states at the Fermi level Ep, and 
it is strongly anisotropic in its crystal and electronic 
structures. [Q-|8j In another important way, it is dis- 
tinct from the other intermetallic superconductors: 
it derives its high superconducting critical tempera- 
ture T c from extremely strong coupling from only a 
small fraction of the phonons to specifically a limited 
part of the Fermi surface (FS). @,| Other 
diborides, at least so far, are disappointing with re- 
gard to their superconductivity - a study of TaB2 
reveals that the presence of Ta d electrons in the 
valence band region results in a completely different 
electronic structure, especially the states at the 
Fermi level, and accounts for its lack of bulk super- 
conductivity. 

Although the properties of MgB2 appear to be 
described consistently, and reasonably accurately, 
by a Fermi liquid picture based on the band struc- 
ture calculated in the local density approximation 
(LDA), there have been few opportunities to make 
detailed quantitative comparison with experimental 
data. Therefore there has been a broad range of al- 
ternative suggestions. Imada has suggested strong 
interband Coulomb exchange processes p^] . Fu- 
rukawa has raised the possible importance of half- 
filled p z bands [fj3[ . Hirsch and Marsiglio Q sug- 
gest that hole-undressing by Coulomb interaction in 
nearly filled bands is responsible. Referring to opti- 
cal data on oriented films, Marsiglio |l5| has sug- 
gested that coupling via a hi gh energy electronic 
mode is plausible. Baskaran |l6| has revived the 



Pauling resonating-valence-bond picture of benzene 
as having possible application to the graphene layers 
of boron in MgB2. Nonadiabatic processes strongly 
affecting the occurrence of superconductivity have 
been put forward, by Alexandrov based on penetra- 
tion depth data p"^ ] and by Cappelluti et al. on 
the basis of several experimental results that are 
not readily understandable in terms of conventional 
(isotropic) Eliashberg theory. |l^] 

Although the crystal structure of MgB2 is quite 
simple, it is strongly layered and electronic charac- 
teristics are predicted to be strongly anisotropic. So 
far few single crystals have been obtained, so single 
crystal optical conductivity data is not available. In 
addition, the strong "anisotropics" mentioned above 
necessitate a two- band model 0,|[^9 20 or a fully 
anisotropic treatment |L0| to account properly for 
the data, and such interpretations are only recently 
being carried out. Photoemission data on single 
crystals allow direct comparison, and the agreement 
between band theory and angle-resolved photoemis- 
sion spectroscopy is good pl|^2[| down to the scale 
of a few tenths of an eV. 

The recently reported observation of de Haas- van 
Alphen (dHvA) oscillations on single crystals by Yel- 
land et al. p3[ | is a crucial development that provides 
the opportunity for detailed evaluation of LDA pre- 
dictions. These measurements detected three fre- 
quencies (extremal Fermi surface areas F) and in- 
formation about the orbit-averaged electron-phonon 
effective mass m* and scattering time r. Compar- 
ing to reports of Elgazzar et al. [plf, who used the 
augmented spherical wave (AS Wj method within 
the spherical potential approximation, Yelland et al. 
concluded that discrepancies with band theory are 
40-80% in the FS areas. This discrepancy seems 
large enough to suggest the occurrence of impor- 
tant correlation effects beyond LDA. If this is true, 
our current understanding of the properties of MgB2 
might need revision. Mazin and Kortus have pre- 
sented orbital areas [B5| that are close enough to the 



observed values to give confidence in the band pic- 
ture. 

In this paper we provide a comparison of care- 
ful calculations of extremal areas and band masses 
for MgB2. To provide comparison and contrast in a 
related diboride, we provide similar information for 
the isostructural refractory diboride ZrB2, which has 
recently been reported [[26j to superconduct at T c 
= 5.5 K whereas earlier searches observed no super- 
conductivity. J27J The extensive dHvA data available 
for ZrB2 makes this a particularly useful system to 
study. We find, with one possible exception that we 
discuss, that LDA predictions seem to provide an ex- 
cellent description for ZrB2. For MgB2 there is some 
disagreement in FS areas that can be accounted for 
by a shift of tt (B p z ) bands with respect to a (B 
sp x py) bands by 240 meV and readjustment of the 
"Fermi energies" of each of these bands by ~ ± 120 
meV. This result is roughly consistent with that of 
Mazin and Kortus. Q Such shifts will lead to quan- 
titative, but most likely not qualitative, corrections 
in the extent explanation of MgB2 superconductiv- 
ity. 



II. CALCULATIONAL METHODS 



suited in agreement with the FLAPW result. On the 
other hand, the standard FPLO-basis set resulted in 
agreement with the WIEN97 code for ZrB2, suggest- 
ing that the charge density for the latter compound 
is less anisotropic than in MgB2. This implies as 
well that the discrepancy of the results of Elgazzar 
et al. [p4f is due to the spherical approximation of 
the potential they used. 

The second item is the dense k point sampling that 
is necessary to obtain the fraction of a band holes ac- 
curately, and hence the charge density and potential, 
a point noted by Mazin and Kortus. [ p5| Increasing 
the number of inequivalent, equally spaced points in 
the irreducible BZ from 2000 (which already would 
be considered to be a fine mesh) to ^16600 results 
in changes in area of the smallest orbits (see Table 
I) by 2-3%. Because of such sensitivities, we quote 
areas and masses only to 2+ significant digits. Un- 
like non-sphericity and k point sampling, the choice 
of exchange-correlation potential makes no physi- 
cal difference. Using the LAPW code, we checked 
the eigenvalues of the unoccupied a states at T and 
A. Measured relative to the Fermi level, these dif- 
fered by no more that 2.5 meV between the LDA 
of Perdew and Wang |K| and the generalized gradi- 
ent approximation of Perdew, Burke, and Ernzerhof. 



The precision and consistency of the calculations 
using different methods is a concern, in light of the 
apparent discrepancy between the band structure re- 
sults of Elgazzar et al. and the data of Yelland et 
al. Therefore we have applied two full potential, 
all-electron methods of calculation that have pro- 
duced equivalent results for several other systems. 
p8| , p9| One method is the full potential linearized 
augmented plane wave (FLAPW) method |3(| as 
implemented in the WIEN97 code. (3l) The other 
method is full potential local orbital code (FLPO) 
p2[ based on local orbitals optimized to minimize 
the electronic energy. 

There are two computational details that require 
attention, especially for MgB2, in order to obtain the 
precise predictions of band theory. Both are related 
to the existence of small volume Fermi surfaces [the 
smallest contains ~3% of the Brillouin zone (BZ) 
volume], for which small shifts of the band edge 
make appreciable differences. The first issue is that 
the non-spherical nature of the charge density (re- 
lated to the distinctly different contributions to the 
density from a and tt states) and potential is impor- 
tant. Comparing the FPLO method with basis func- 
tions optimized as usual (to minimize the energy) 
with the FLAPW results revealed that the k z disper- 
sion of the a band along T-A was slightly different 
for the two methods. The reason is that the exten- 
sion parameter Xq' 1 for the FPLO basis functions is 
optimized only with respect to the main quantum 
number n and the angular momentum I at a given 
site. Increasing the flexibility of the FPLO basis re- 



III. RESULTS AND COMPARISON WITH 
EXPERIMENT 

A. ZrB 2 

As in MgB2 (but for a different reason) it is not 
easy to prepare single crystals of ZrB 2 due to its 
melting temperature above 3300 K. However, single 
crystal studies have been reported. dHvA data have 
been provided and analyzed for ZrB2 by Tanaka and 
Ishikawa, |35|,[36| and the data are similar to those re- 
ported for the isostructural and isovalent sister com- 
pounds TiB 2 HQ and HfB 2 . §| 

The band structure has been presented previ- 
ously, based on a variety of computational meth- 
ods, [^9|-^3| each involving simplifications that we 
avoid. Our calculated band structure and density of 
states (DOS), calculated for a =3.170 A, c =3.532 A, 
are similar to those presented earlier, so we do not 
present them here. The Fermi energy lies in a "pseu- 
dogap," which is due to relatively large velocities in 
the region of E^ rather than any semimetallic over- 
lapping of bands (where there would be small FSs 
and low velocities). The occupied B ir states are 
mostly in the -4 eV to -2 eV region (E^ is taken 
as the zero of energy), while a states are spread 
throughout the valence band region. Due to strong 
hybridization of the Zr 3d states with the B 2p states, 
the cr-7r distinction is not as clear in ZrB 2 as it is in 
MgB2. The Zr 3c? DOS shows that considerably less 



than 40% of the 3d DOS is occupied (which would be 
the case four of ten possible 3d electrons) suggesting 
Zr— >B charge transfer. 
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FIG. 1. Total and projected density of states for 
ZrB2. The projections show that there is no distinctive 
character of states at Ef, rather they are a combination 
of Zr Ad and B n and a. 

The Fermi surfaces, shown in Fig. g, consist of a 
K-centered barbed ring 1Z with threefold symmetry, 
and an A-centered dumbbell V with sixfold symme- 
try. We use the orbit designations of Tanaka and 
coworkers. |3^|3(j] For the field along (0001), the 1Z 
surface gives a rounded triangular orbit (v) encir- 
cled by a nearly circular orbit (£), while the T> sur- 
face gives a circular orbit (e) around its waist and 
a smoothed-star orbit (fi) at each end. These cross 



4 In addition, we 
ong (1010), which 



sections are also shown in Fig. 
consider the f3 orbit for field at 
is the cross section of the barbed ring V in Fig. |2j. 
These Fermi surfaces, which are broadly consistent 
with those used by Tanaka and coworkers to inter- 
pret their dHvA data. The Fermi surfaces presented 
by Shein and Ivanovskii (Fig. 1 in Ref. p^]) are 
quite different. The origin of this difference is un- 
clear, since the full potential linear muffin-tin orbital 
method they use should give the same results as our 
methods. We note that their value N(E^) = 0.163 
cV _1 differs considerably from our results and most 
of the previous calculations. |3^ , ^0| , ^6| 




FIG. 2. Calculated Fermi surfaces (left) and selected 
cross sections (right), for ZrB2. The labels are provided 
as used in the text. 

The general shape of the density of states, and the 
low value of N(E^ ), was confirmed by xray photoc- 
mission spectroscopy measurements. |39|l The most 
noteworthy feature, in light of a new report p(| of 
T c = 5.5 K, is the small calculated value of N(E^) 
= 0.26 /eV-cell, corresponding to a bare linear spe- 
cific heat coefficient jb — 0.61 mJ/mole-K 2 . [The 
value of N(Ei?) is somewhat sensitive to the qual- 
ity of the calculation.] The reported experimental 
Q value of 7 = 0.47 mJ/mole-K 2 , which should in- 
clude electron-phonon enhancement, is smaller than 
our bare band structure value. 

The low value of N(E,p) itself and the weak cou- 
pling |l5| seem inconsistent with T c = 5.5 K; it 
would require very large electron-phonon matrix el- 
ements which has been shown not to be the case in 
another transition metal diboride TaB2. (T^| Since 
all information is consistent with the observed su- 
perconductivity arising from a minority phase in the 
sample, and the calculations seem to make make su- 
perconductivity unlikely, we conclude that ZrB2 it- 
self is not superconducting. These facts, and results 
presented below, suggest that remeasurement of the 
heat capacity on additional samples may be called 
for. In any case, the occurrence (or not) of super- 
conductivity in ZrB2 is peripheral to the intent of 
this paper. 
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FIG. 3. The areas F(E) of selected extremal constant 
energy surfaces in ZrB2 for field along (0001), versus en- 
ergy (Ef = 0). Sloping lines show our calculated values, 
experimental numbers are denoted by horizontal lines, 
and the circles denote the energy where perfect agree- 
ment occurs. The three smaller are in good agreement 
with experiment; the discrepancy for the large surface is 
discussed in the text. 

The calculated and observed areas are given in Ta- 
ble I. Percentage discrepancy or even absolute dif- 
ferences do not give the most physical indication of 
the level of (dis)agreement, especially when areas get 
small (see MgB2, below). To indicate more clearly 
the implication of the discrepancy between theory 
and observation, in Fig. || the orbit areas F(E) ver- 
sus energy are provided. The observed areas are 
shown as horizontal lines, which allows one to read 
off the shift in any given band to bring its area into 
agreement with the observed area. The band masses, 
proportional to the derivative dF/dE, are also read- 
ily obtained from such curves. 

For the orbits e,n,v in the range 1.8 kT < F < 
3 kT, the agreement is excellent, probably within 
the total numerical precision. For the £ orbit, the 
agreement is not good: the reported area is 10% 
smaller than we calculate (this amounts to 1.4% of 
the BZ basal plane area). Moreover, it is unlikely 
that even this relatively modest disagreement can 
be achieved by a shift of the band, because such a 
shift destroys agreement for the v orbit and results 
in a change of topology of the constant energy sur- 
faces. To investigate this discrepancy further, we 



have calculated the (3 orbit around the 1Z surface, 
which passes through both the £ and v orbits. This 
area is in excellent agreement with experiment (see 
Table I). It may be possible, but seems unlikely, that 
this disagreement could be reconciled by some non- 
rigid shift(s) of band(s). It is possible that this dis- 
crepancy is due to experimental conditions: the £ 
orbit is the largest, giving the fewest oscillations to 
ht to an oscillatory form, hence larger uncertainty 
in the result. An experimental reinvestigation of the 
dHvA frequencies is under way to either confirm 
or resolve this discrepancy. 

B. MgB 2 

To facilitate understanding of notation, we iden- 
tify the orbits by their B character (a or tt) and by 
the point in the BZ around which they are centered. 
The a Fermi surfaces, pictured in several previous 
publications, [||](| are two concentric fluted cylinders 
oriented along T-A, which give rise to extremal or- 
bits Cp, <7p, o\, a\ for magnetic field along (0001). 
(S, L denotes small, large.) The it bands give rise to 
7tr and tta for field along (0001), ttm for field along 
(1100), and tt l for field along (1000). 

The calculated areas and band masses are given in 
Table II, with comparison to the three orbits of Yel- 
land et al. |23f| assuming the same correspondence 
of observed orbit and calculated orbit. For refer- 
ence, an area of 1 kT corresponds to 2% of the area 
Abz of the basal plane of the BZ. For the three 
observed orbits, the calculated areas are 0.30±0.04 
kT larger than observed, i.e. a discrepancy equal 
to 0.6% of Abz- Our calculated areas are in good 
agreement (usually close to significant digits) with 
those of Mazin and Kortus. |2^] There are differ- 
ences compared to the areas presented by Elgazzar 
et al. Q| (most of our areas are ~20% smaller), pre- 
sumably due to the approximations made in their 
augmented spherical wave method. 

In Fig. |J the orbit areas F(E) versus energy are 
provided. The band masses mi, = (H 2 /2Tr)dF/dE 
are seen to be insensitive to the position of the Fermi 
energy except for the two largest n orbits. These 
values are in good, but not quite perfect, agreement 
with those of Mazin and Kortus. The observed ar- 
eas (Table II) are shown as horizontal lines, which 
allows one to read off the shift in any given band 
to bring its area into agreement with the observed 
area. (Due to the difference in the a and 7r densities 
of states at Ef , charge balance requires a readjust- 
ment of Ef by ^20 meV.) The required energy shift 
is -115 meV for both cr§ and cr^, and +125 meV for 
7T£,, both of which reduce the sizes of the hole and 
electron Fermi surfaces. These required shifts imply: 
(1) there is a relative shift of the a bands relative to 
the 7T bands by ~ 240 meV presumably due to cor- 
relation effects beyond LDA, and (2) the k z disper- 



sion of the a bands along the T-A line is described 
correctly in band theory (the same a band shift is 
required at T and at A). The data give the areas for 
only a single tube (the smaller one), but since the 
two er are degenerate along T-A, both bands must 
shift together. 
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FIG. 4. The areas F(E) of selected extremal constant 
energy surfaces as in Fig. H, but for MgB2 . Bands shifts 
of about -115 meV for the a bands and +125 meV for the 
7r bands bring the calculated areas into perfect agreement 
with the data. 

The observed effective masses include an enhance- 
ment due (primarily) to electron-phonon coupling, 
which is given in Table II and is obtained from 
to* = mt,(l + A). For both a orbits the enhance- 
ment thus derived is X a =1.3; for the tt orbit it 
is A 7r =0.5. These a orbit enhancements are notice- 
ably larger than the average values over each surface 
A (j = 0.9 ± 0.1, K = 0.4 ± 0.1 obtained from so- 
lution to the anisotropic Eliashberg equations with 
first principles band theory input by Choi et al. |h]] 
The difference particularly for the a bands might 
indicate there is some feature of strong coupling in 
MgB2 yet to be understood, or possibly that the 
a band shift necessary to provide agreement with 
orbit areas, when taken into account in the calcu- 
lations of electron-phonon coupling, will ameliorate 
this discrepancy. The disagreement is small enough 
that the general picture - the quasiparticlc band 
structure and electron-phonon coupling determined 
within conventional (LDA) band theory, then T c de- 



termined by Eliashberg theory - seems well justified. 



IV. DISCUSSION 

For ZrB2 the only discrepancy between calculated 
and observed areas occurs for the largest orbit that 
was observed, the overall agreement between the cal- 
culated and the observed frequencies is of the same 
order as the achievable numerical accuracy. Because 
of the strong hybridization between the Zr 3d and 
the B 2p states, ZrB2 has a much more isotropic elec- 
tronic structure than does MgB 2 . As a result, any 
renormalization and shift of bands with respect to 
each other seems to be negligible for this compound, 
and the LDA single particle picture description is 
well justified. 

While the observed dHvA areas of MgB2 are read- 
ily understood in terms of the calculated Fermi sur- 
faces and the inferred electron-phonon coupling con- 
stants are reasonable, complete agreement of the FS 
areas with the data requires a shift of the a bands 
by 115 meV downward, and a shift of the tt bands 
upward by 125 meV. The necessary shift results in a 
band edges at T and A that are reduced 0.38 eV — > 
0.25 eV, 0.76 eV -» 0.63 eV, respectively. The vol- 
ume of the a s tube, proportional to the average of 
the two areas in Table II, then is about 30% smaller 
than given by the band calculation, so the number of 
holes decreases proportionally (from 0.146 to 0.106). 

We can suggest at least two possible causes of 
this "beyond LDA" correction to the band structure. 
One possibility is related to the observation that the 
charge in the a bands is confined to the two dimen- 
sional B sheet of graphene structure, nearly filling 
those bonding states. Being of more localized char- 
acter than the tt bands, the exchange potential in 
LDA may be less accurate than for the tt electrons; 
a better exchange potential would be larger in mag- 
nitude (and attractive), lowering the a bands with 
respect to the tt bands. A related viewpoint of the 
same physics is that there is a larger (spurious) self- 
interaction for the a states than for the tt states in 
LDA. 

Another possible correction could arise from resid- 
ual hole-hole Coulomb interactions in the cr-hole gas. 
There is some analogy with a related situation in 
metallic, ferromagnetic Ni, where there are ~0.6 3d 
holes/Ni in the minority bands. Unlike in MgB2, 
in Ni there is presumably a relatively large value of 
the ratio U/W (intra-atomic Coulomb interaction 
U, 3d bandwidth W) making Hubbard-type correla- 
tions of some relevance. MgB2 has broad a bands, 
so Hubbard-like correlations should not be the prob- 
lem. The strong 2D character of the a band holes 
may enhance many body corrections. It will require 
further work to determine whether it is one of these 
mechanisms, or perhaps some other, that is respon- 
sible for the band structure corrections. 



V. SUMMARY 

The rather close correspondence of the LDA Fermi 
surfaces of MgB2 supports the prevalent picture of 
superconductivity that is based on the LDA start- 
ing point: very strong electron-phonon coupling of 
a band hole states to certain phonons (the E2 9 
branch). Of the suggested alternatives cited in the 
Introduction, the comparison we have made tends 
to rule out several of them. Although the a band 
Fermi energy requires some adjustment to account 
for the observed dHvA orbital areas, it remains 
around 0.5 eV, large enough to render non-adiabatic 
processes unimportant in superconducting pairing. 
The half-filled tt ideas and resonating valence bond 
state idea certainly get no support from the corre- 
spondence of dHvA data to the band theory results. 
Other experimental data (penetration depth, spe- 
cific heat, isotope shift, superconducting gap and T c 
itself) are being shown in many studies to require an 
anisotropic (two band) model, and that most data 
seem to be consistent with such models based in de- 
tail on band theory results. Single crystal data, of 
which the dHvA data Q are some of the first, will 
serve to clarify these issues further. 
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TABLE I. Calculated de Haas-van Alphen areas and 
masses of ZrB2, compared to the experimental data of 
Tanaka and Ishizawa. Areas are quoted in kTesla. 



Orbit 


Fexp 


FcaZc 


m b 


£(0001) 


1.81 


1.84 


-0.38 


At(0001) 


2.46 


2.43 


-0.60 


i/(0001) 


2.84 


2.92 


-0.41 


£(0001) 


6.09 


6.78 


0.51 


/3(1010) 


0.300 


0.305 


0.103 



TABLE II. Calculated de Haas-van Alphen parame- 
ters of MgB2 compared to the experimental data. Ar- 
eas are quoted in kT. Column 1: the calculated FPLO 
areas, using the LDA potential calculated using 16221 
k points in the irreducible BZ; values in parentheses are 
from FLAPW using the CCA exchange-correlation. Col- 
umn 2 and column 4: data from Yelland et al. Column 
3: band mass, in units of electron mass (negativ masses 
stand for holes). Column 5: orbit-averaged mass en- 
hancement A = |ra*|/|ra(,| — 1. 



Orbit 






mb 


\m* 


A 
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0.78 (0.79) 


0.54 


-0.25 


0.57 


1.3 




1.65 (1.67) 




-0.57 






tit 


34.5 




1.87 






*% 


1.83 (1.81) 


1.53 


-0.31 


0.70 


1.3 




3.45 (3.46) 




-0.64 








30.6 




-0.93 






km 


0.45 




-0.25 






7VL 


3.03 


2.69 


0.32 


0.47 


0.5 



